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ABSTRACT 


Shock  wave  thicknesses  and  density  ratios  have  been  measured  in  helium, 
argon  and  nitrogen  by  means  of  the  electron  beam  fluorescence  method,  over  the 
range  1.5  <  M  <  17.4,  in  a  low  density  wind  tunnel.  The  sh ;cl;  thicknesses  in 
argon  and  helium  agreed  well  with  Mott-Smith  theory  at  the  higher  Mach  numbers, 
and  were  between  Navier-Stokes  and  Mott-Smith  theory  at  the  lowest  Mach  number. 
In  nitrogen,  the  measured  shock  thicknesses  were  considerably  greater  than  the 
predictions  of  Navier-Stokes  theory. 

Measured  density  ratios  across  the  shock  wave  were  in  good  agreement  witn 
theory,  in  the  lower  density  flows.  Poorer  agreement  was  found  at  higher  flow 
densities,  leading  to  estimates  of  upper  bounds  for  the  range  of  linear  varia¬ 
tion  of  fluorescence  intensity  with  gas  density.  Measured  density  ratios  for 
shocks  produced  in  divergent  free-jet  flows  were  found  to  be  in  better  agree¬ 
ment  with  theory  after  a  viscous  curvature  correction  was  applied. 


1.0  INTRODUCTION 


Several  different  techniques  have  been  employed  successfully  to  measure 
the  thickness  and  internal  structure  of  shock  waves.  The  earliest  determina¬ 
tions  of  shock  wave  thickness  were  r  _e  by  means  of  the  optical  reflectivity 
technique  (see  Linzer  and  Hornig,^  and  the  references  cited  therein).  At 
about  the  same  time  the  optical  reflectivity  technique  was  being  perfected, 

measurements  of  the  internal  structure  of  shock  waves  were  carried  out  by 

2  3 

means  of  free-molecule  probes.  5  Results  obtained  by  these  two  techniques 

and  comparisons  with  theory  are  discussed  in  the  survey  papers  by  Sherman  and 
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Talbot  and  Talbot. 

mcc-  r2ccr.t  1  ■» cKa/'I-  uauo  c frnr ♦‘mta  ti£ivG  been  by 

electron  beam  techniques.  The  idea  of  using  an  electron  beam  dates  back  at 
least  to  the  work  of  Ballard  and  Venable,^  who  observed  the  change  in  trans¬ 
mission  of  an  electron  beam  directed  across  a  shock  tube,  due  to  the  passage 
of  a  shock  wave.  Only  a  few  data  points  were  obtained,  however,  and  these 
were  subject  to  rather  large  experimental  uncertainty.  Recently  Russell^  has 
succeeded  in  improving  the  accuracy  of  the  method  considerably,  mainly  by  the 
use  of  a  large-diameter  low-density  shock  tube,  and  has  reported  new  data  on 
the  density  thickness  of  argon  shock  waves.  Similar  work  has  been  reported 

g 

by  Schultz-Grunow  and  Frohn. 

Another  method  of  using  an  electron  beam  for  probing  shock  waves  has  been 
9 

developed  by  Camac,  who  observes  the  electrons  scattered  out  of  the  beam  by 
collisions,  rather  than  the  portion  of  the  beam  which  is  transmitted  unatten¬ 
uated  across  the  shock  tube.  With  this  method,  Camac  has  obtained  shock 
thickness  data  for  argon,  nitrogen  and  carbon  dioxide.  His  argon  results  are 
in  reasonably  good  agreement  with  those  of  Russell. 

Still  another  technique  involving  the  use  of  an  electron  beam  is  available 
for  studying  shock  wave  structure.  This  is  the  electron  beam  fluorescence 


3. 

fairly  comprehensive  picture  of  the  problem,  especially  for  monatomic  gases- 
Our  main  purpose  in  carrying  out  the  present  series  of  experiments  has  in  fact 
been  to  use  the  shock  wave  flow  as  a  testing  ground  for  assessing  the  accuracy 
and  limitations  of  the  electron-beam-fluorescence  technique,  although  of  course 
in  the  process  of  so  doing  we  have  attempted  to  extend  the  available  data  and 
obtain  certain  new  results. 

2,0  DESCRIPTION  OF  THE  METHOD 
2.1  Equipment 

The  theory  underlying  the  electron-beam  fluorescence  technique  is 
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described  in  the  papers  by  Muntz  and  Muntz  and  Marsden,  We  shall  not  go 
into  an  extended  description  here,  and  we  refer  the  reader  to  these  papers  for 
details. 

The  experimental  arrangement  used  in  the  present  tests  is  shown  in 

Fig,  1.  The  electron  beam  source  is  a  commercial  Alloyd  electron-beam  welder, 

rated  at  30  kV  and  150  mA.  (In  the  tests,  beam  currents  in  the  range  1,5  - 

10  mA  were  used,  at  30  kV„)  Downstream  of  the  exit  orifice  of  the  beam  chamber 

a  magnetic  X-Y  deflection  coil  system,  also  manufactured  by  Alloyd,  was  used  to 

direct  the  beam  into  the  desired  trajectory.  The  wind  tunnel  static  pressure 

was  typically  in  the  10  ^  Torr  range,  and  a  3-inch  oil  diffusion  pump  attached 

to  the  beam  chamber  was  adequate  to  keep  the  beam  chamber  pressure  below 
-L 

2  X  10  Torr  when  a  0,040  inch  diameter  exit  orifice  was  used. 

Shock  waves  were  produced  in  the  flow  by  means  of  a  shock-holder. 

A  water-cooled  "choker"  device  downstream  of  the  shock-holder  was  positioned 
so  that  the  shock  wave  stood  slightly  ahead  of  the  shock-holder.  Slots  about 
1/4  inch  long  were  cut  in  the  sides  of  the  shock-holder,  so  that  the  flow  could 
be  observed  some  distance  downstream  of  the  entrance.  The  shock  waves  could  be 
observed  visually  in  the  darkened  tunnel,  because  the  primary  electron  beam. 
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together  with  the  "halo"  around  it  (presumably  produced  by  secondary  electrons) , 
illuminated  a  substantial  area  of  the  shock  wave.  All  of  the  shock  waves  studied 
appeared  to  be  quite  plane,  although  no  precise  measurements  of  shock  curvature 
were  attempted.  (Some  discussion  of  shock  curvature  effects  is  given  later  in 
this  report.)  The  shock-holder  choker  had  a  0.080  inch  diameter  collimating 
holes  through  which  the  electron  beam  passed;,  at  a  slight  angle  to  the  axis  of 
the  flow,  and  upstream  of  the  shock  wave  the  beam  was  collected  by  a  water- 
cooled  cup.  The  current  to  this  cup  was  used  in  a  feedback  loop  to  regulate 
the  voltage  on  the  bias  grid  of  the  electron  gun,  in  order  to  maintain  a  con¬ 
stant  beam  current. 

The  optical  detection  system  is  also  illustrated  schematically  in 
Fig.  1,  The  entrance  slit  of  a  Jarrell-Ash  3.4  meter  Ebert  spectrograph  was 
imaged  at  approximately  unity  magnification  perpendicular  to  the  electron  beam, 
by  a  48  inch  focal  length  spherical  mirror  and  a  plane  mirror,  as  shown.  The 
slit  image  was  scanned  along  the  electron  beam  by  rotating  the  plane  mirror, 
and  the  output  of  a  10-turn  potentiometer  linked  to  the  motor-driven  micrometer 
drive  of  the  plane  mirror  was  coupled  to  the  X-axis  of  an  X-Y  recorder.  The 
radiation  from  the  beam  was  detected  by  an  uncooled  IP-21  photomultiplier  tube, 
whose  output  was  amplified  by  a  Model  600A  Keithley  electrometer  and  fed  to  the 
Y-axis  of  the  recorder.  A  stack  of  ten  67.5-Volt  batteries  was  used  to  power 
the  photomultiplier.  While  batteries  probably  have  no  particular  advantage  over 
a  regulated  power  supply  used  with  a  voltage  divider,  they  have  proven  to  be 
stable,  trouble  free,  and  inexpensive.  In  almost  all  cases,  the  photomultiplier 
was  placed  directly  behind  the  50  micron  entrance  slit,  so  that  the  spectrometer 
was  not  used  as  such.  The  only  spectral  resolution  of  the  radiation  in  the 
actual  shock  wave  measurements  was  that  afforded  by  a  Corning  7-59  blue  filter 
placed  ahead  of  the  entrance  slit  in  the  argon  and  nitrogen  measurements.  The 
filter  was  not  used  for  the  helium  measurements.  However,  spectral  analyses 
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of  the  radiation  were  carried  out  under  both  static  and  flow  conditions.  The 
principal  source  of  the  observed  radiation  in  nitrogen  was  found  to  be  the 

O 

3914  A*  (0*0)  band  of  the.  first  negative  system  of  •  In  argon*  the  radiation 

o 

comes  principally  from  lines  of  AI  and  ArT  between  3000  and  4000  A,  although 

there  are  a  number  of  lines  in  the  red  which  were  excluded  by  the  filter.  In 
3  3 

helium  the  5015  A  3  P  and  4471  A  4  D  lines  are  responsible  for  most  of  the 
radiation;,  the  latter  line  having  about  5%  the  intensity  of  the  former. 

2.2  Spatial  Resolution 

The  resolution  of  the  optical  system  was  estimated  to  be  about 
0.002  inch,  measured  along  the  beam  axis,  when  a  0.050  mm  entrance  slit  was 
used.  This  estimate  was  arrived  at  by  scanning  the  slit  image  across  a 
knife-edge  located  within  the  wind  tunnel.  The  knife-edge  was  illuminated 
from  behind,  and  the  resolution  was  taken  arbitrarily  to  be  the  distance  of 
slit-image  travel  corresponding  to  the  increase  from  10%  to  90%  of  the  total 
light  intensity  observed  by  the  photomultiplier  as  the  slit-image  traversed 
the  occulting  knife-edge. 

An  additional  factor  contributing  to  the  overall  resolution  is  the 

convection  downstream  of  the  radiating  atoms  and  molecules.  Some  measurements 

of  this  effect  were  carried  out  by  passing  the  electron  beam  vertically  upward 

through  a  0.018  inch  scmi-eircular  aperture  (formed  by  half-plugging  a  0.P40 

inch  hole,  with  the  flat  edge  of  the  aperture  upstream)  in  a  water-cooled 

wedge  whose  upper  surface  was  aligned  with  the  flow  issuing  from  a  Mach  4 

nozzle.  In  the  free  stream  region  above  the  wedge  shock  wave,  the  beam 

fluorescence  was  scanned  with  and  without  flow,  in  order  to  determine  the 

convective  displacement  of  the  radiation  profile.  This  displacement  x  is 

approximately  related  to  the  lifetime  x  of  the  excited  state  through 

x  =  Ux,  where  U  is  the  stream  velocity.  On  the  basis  of  the  measured 

-8  + 

lifetime  T  =  6.5  X  10  sec  of  the  (0,0)  band  of  Ng  (see  Ref.  15)  the 
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displacement  x  in  ^  should  have  been  about  0.0015  inch,  whereas  we  observed 

x  <  0.002  inch,  with  an  experimental  uncertainty  of  about  0,002  inch.  No  effort 

was  made  to  determine  the  value  more  precisely.  Similar  small  displacements 

were  observed  for  the  total  radiation  in  argon  (with  blue  filter),  and  the 

«  1 

spectrally  resolved  5015  A  3~P  helium  line.  On  the  other  hand,  the  spectrally- 

o  3 

resolved  radiation  from  the  4471  A  4  D  helium  line  showed  a  significant 

(x  »  0.006  inch)  downstream  shift. 

The  observed  radiation  profiles  are  shown  in  Figs.  2a  and  2b.  All 

represent  data  taken  under  flow  conditions.  It  will  be  noted  that  the  profiles 

obtained  from  argon,  nitrogen  and  the  3^P  helium  line  are  quite  symmetrical. 

3 

The  4  D  helium  line,  on  the  other  hand,  shows  a  pronounced  asymmetry  in  the 
downstream  direction.  Of  particular  interest  in  these  profiles  is  the  inten¬ 
sity  of  the  radiation  in  the  "halo"  surrounding  the  primary  beam.  This  halo 
is  believed  to  be  due  to  excitation  of  the  gas  by  secondary  electrons  resulting 
from  ionizing  collisions  in  the  main  beam,  rather  than  to  excitations  by  primary 
electrons  which  have  undergone  small-angle  scattering  out  of  the  main  beam. 

These  secondary  electrons  will  be  of  predominantly  low  energy,  and  in  fact  many 
of  them  may  have  insufficient  energy  to  excite  the  observed  N^+  radiation  (see 
Muntz  and  Marsden,  and  Mott  and  Massey^  ),  Thus  they  will  have  a  short  mean 
free  path,  and  would  be  expected  to  leave  the  primary  beam  by  diffusive  processes. 
The  measurements  of  halo  intensity  shown  seem  to  bear  out  these  predictions.  The 
magnitude  of  the  halo  in  different  gases  would  depend  on  the  particular  cross 
sections  for  the  various  important  processes.  The  smallest  relative  nalo  inten¬ 
sity  occurs  in  helium,  with  nitrogen  only  a  little  larger,  while  argon  has  a 
quite  large  halo  intensity.  Some  measurements  were  made  at  different  gas  den¬ 
sities,  covering  about  a  factor  of  3,  and  while  in  general  the  halo  seemed  to 

increase  in  relative  intensity  with  gas  density,  the  effect  was  not  pronounced, 

3  1 

The  4  D  line  of  helium  shows  somewhat  more  halo  than  the  3  P  line;  however,  the 
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4  D  line  is  only  about  5%  the  intensity  of  the  3  P.  The  measured  shifts  of 

these  two  lines  of  helium,  and  the  halo  intensity,  are  not  understood,  as  they 

17 

do  not  seem  to  correlate  with  the  theoretical  lifetimes  which  have  also  been 

15  o 

verified  by  experiment .  It  must  be  remembered  that  a  537  A  resonance  line 

also  originates  from  the  3*P  level,  and  at  the  density  and  temperatures  employed 
here  this  radiation  will  be  strongly  imprisoned,  so  that  Cue  Iiretime  of  the 
31?  level  should  oe  given  by  decay  through  the  5015  A  radiation  only. 

The  halo  radiation  is  of  concern  in  the  shock  wave  measurements  because 
one  must  look  through  the  halo  when  observing  the  primary  beam,  even  if  the 
spectrograph  slit  height  is  chosen  no  larger  than  the  primary  beam  diameter. 

In  point  of  fact,  it  is  preferable  for  the  shock  wave  measurements  to  choose 
the  slit  height  somewhat  greater  than  the  primary  beam  diameter,  because  in  this 
way  a  small  deviation  from  parallelism  between  the  beam  and  the  path  followed 
by  the  slit  image  will  not  produce  large  errors;  this  procedure,  however, 
increases  the  amount  of  halo  radiation  included  with  that  of  the  main  beam. 

Now,  if  the  halo  radiation  intensity  varies  linearly  with  density,  in  the  same 
way  as  does  the  primary  beam,  it  should  not  present  much  of  a  problem  in  making 
shock  wave  density  profile  measurements.  It  appears  from  our  data  that  the 
variation  is  approximately  linear,  but  the  evidence  is  not  conclusive.  But 
even  if  the  halo  intensity  does  vary  linearly  with  density,  it  still  appears 
to  present  a  problem  in  making  rotational  temperature  measurements  in  nitrogen. 
This  is  discussed  in  the  companion  reports. ^ 

The  combination  of  uncertainties  due  to  the  resolution  of  the  optics 
and  the  downstream  convection  of  excited  species  leads  to  the  estimate  of 
approximately  0.003  inch  probable  error  in  the  location  of  a  point  in  the  flow. 
Since  the  shock  waves  reported  here  have  maximum  slope  thicknesses  greater  than 
0.030  inch,  there  should  be  no  measurable  distortion  of  either  the  shock  wave 
profile  or  the  maximum  slope  thickness,  as  pointed  out  by  Russell. ^ 
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2.3  Free  Stream  Flows 

Both  contoured  nozzles  and  free  jet  expansions'5^  were  used  to  pro¬ 
duce  the  test  flows.  The  gases  were  unheated,  so  that  all  stagnation  temper¬ 
atures  were  room  temperature  (~  297°K).  All  flows  were  calibrated  by  measuring 
stagnation  and  impact  pressures8  and  assuming  isentropic  flow.  The  impact 
pressures  were  corrected  experimentally  for  viscous  effects  by  testing  a  range 
of  sizes  of  open-ended  and  source-shaped  probes,,  and  extrapolating  to  infinite 
probe  diameter.  In  the  case  of  the  free  jet  flows,,  a  correction  was  also  made 
for  shock-standoff  distance.  This  correction  is  necessary  because  an  axial 
gradient  is  present  in  the  free  jet  flow,  and  the  stream  Mach  number  associated 
with  the  probe  position  is  slightly  higher  than  that  of  the  probe  shock  wave. 

Variation  of  the  Mach  number  in  the  free  jet  tests  was  achieved  by 
keeping  the  shock-holder  at  a  fixed  position  (about  6  inches  downstream  of  the 
orifice  plates)  and  changing  the  diameter  of  the  sonic  orifice.  This  procedure 
was  more  convenient  than  varying  the  distance  between  the  shockholder  and  a 
fixed-diameter  orifice.  Orifice  diameters  of  l/4,  l/2,  1  and  4  inches  were 
used.  The  three  smallest  orifice  plates  were  machined  from  l/4  inch  thick 
aluminum  plate,  and  the  orifices  had  1/4  inch  radius  fillets  to  provide  a 
smooth  transition  from  the  upstream  surface  of  the  plate  to  the  exit  of  the 
orifice.  The  4-inch  orifice  was  cut  into  a  0.125"plates  and  the  edges  were 
left  sharp.  At  each  flow  Mach  number,,  for  both  tne  orifices  and  the  nozzles, 
the  stagnation  pressure  pQ  and  hence  the  free  stream  density  was  varied  by 
a  factor  of  4  or  more.  This  variation  produced  negligible  changes  in  M^ 
for  the  orifice  flows,  but  it  did  change  the  exit  Mach  numbers  in  the  nozzle 
flows  because  of  the  change  in  nozzle  boundary  layer  thickness. 

The  No.  4,  Mach  3  nozzle  is  an  axisymmetric  contoured  nozzle  which 
was  fabricated  a  number  of  years  ago.  Although  it  was  not  one  of  the  more 
successful  nozzles  made  for  the  low  density  tunnel,  in  that  it  has  some 
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pronounced  flow  non-uniformities,  it  was  used  in  the  present  tests  because  the 
data  taken  near  Mach  3  using  the  4  inch  sonic  orifice  was  not  felt  to  be  as 
reliable  as  the  other  measurements.  It  was  thought  that  there  may  be  some 
viscous  effects  in  the  flow  produced  with  this  large  orifice.  Unfortunately, 
the  lowest  densities  at  which  the  Mach  3  nozzle  could  be  operated  were  quite 
high,  compared  to  the  rest  of  the  flows,  and  it  appeared  that  the  fluorescence 
technique  develops  some  non-linearities  at  these  densities.  The  No.  11,  Mach 
1.5  nozzle  is  the  monatomic  gas  nozzle  designed  by  Talbot  and  Sherman  for  their 
earlier  shock  wave  studies,  and  it  produces  quite  uniform  flow.  The  data 
obtained  in  the  Mach  3  nozzle  are  considered  to  be  subject  to  more  uncertainty 
than  that  obtained  in  the  Mach  1.5  nozzle. 

Table  I  gives  a  list  of  flow  conditions  representative  of  those  used 

in  the  present  tests.  Also  listed  in  the  Table  are  values  of  the  length 

L  =  n(T  )/pu  which  for  reasons  discussed  later  was  chosen  as  the  normalizing 

length  for  the  measured  shock  wave  thicknesses,  pu  is  the  mass  flux  per  unit 

area  in  the  stream,  which  is  constant  through  a  normal  shock  wave,  and  |i(T  ) 

is  the  viscosity  evaluated  at  the  sonic  temperature  T  =  2Tq/(7+1).  Of  course 
* 

T  is  not  the  temperature  which  actually  exists  at  the  sonic  point  within  a 
shock  wave,  because  the  flow  in  the  interior  of  a  shock  wave  is  generally  not 
locally  adiabatic,  but  it  is  nevertheless  a  useful  reference  temperature. 


3.0  RESULTS 

3-1  Density  Calibrations 

Preliminary  to  the  shock  wave  measurements,  some  static  (no-flow) 
measurements  were  made  in  the  wind  tunnel  of  the  light  intensity  emitted  from 
the  beam  as  a  function  of  gas  density.  It  is  of  course  necessary  to  make  such 
measurements  at  constant  beam  current,  and  it  was  found  that  an  acceptable 
measurement  of  beam  current  could  be  made  only  if  the  beam  was  directed  into 
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a  shielded  cup.  This  cup  consisted  of  a  closed  water-cooled  outer  copper  cyl¬ 
inder,  within  which  was  a  slightly  smaller  copper  cylinder,  electrically 
isolated  but  conduction-cooled  by  boron  nitride  spacers  The  beam  was 
accepted  through  an  aperture  at  one  end  of  the  cup  consisting  of  a  0.102  inch 
diameter  hole  in  the  outer  cylinder*  followed  by  a  0  110  inch  hole  in  the 
inner  cylinder.  The  outer  cylinder  was  held  at  ground  potential,  while  the 
inner  cylinder  was  connected  to  the  current  metering  circuit.  The  object 
of  this  arrangement  was  to  insure  that  the  measured  current  was  due  only  to 
the  primary  high  energy  electrons. 

The  light  intensity  emitted  from  the  beam,  at  a  point  about  1/2  inch 
ahead  of  the  beam  cup,  was  measured  as  a  function  of  gas  density.,  with  the 
beam  current  held  constant.  The  data  indicated  a  reasonably  linear  variation 
of  light  intensity  with  gas  density,  although  the  density  range  which  could 
be  covered  was  not  very  large,  (In  order  to  avoid  spurious  signals  from  the 
afterglow  produced  by  long-lived  metastables,  a  small  flow  was  maintained  by 
the  mechanical  roughing  pump  on  the  tunnel,  and  the  range  over  which  this 
pump  could  operate  fixed  the  density  range  of  the  calibration,)  However,  an 
accurate  calibration  was  not  required  for  the  shock  wave  measurements,  be¬ 
cause  the  density  ratio  across  a  shock  wave  is  known,  and  the  ability  of  the 
electron  beam  technique  to  measure  this  ratio  accurately  at  different  ambient 
densities  is  in  fact  a  quite  satisfactory  test  of  the  linearity  of  the  meas¬ 
urement. 

An  interesting  by-product  of  the  static  calibrations  was  the  obser¬ 
vation  that  although  there  was  considerable  radiation  in  the  beam  halo,  the 
current  in  the  halo  was  essentially  zero.,  This  observation  was  made  by 
traversing  the  entrance  orifice  of  the  shielded  cup  across  the  halo  region. 

The  absence  of  current  in  the  halo  region  is  quite  conclusive  evidence  for 
the  interpretation  that  the  halo  is  produced  by  secondary/  electrons  rather 
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than  by  primary  electrons  which  have  undergone  small-angle  scattering  out  of  the 
main  beam 


32  Shock  Wave  Thicknesses 

Figures  3a  -  3d  show  some  typical  shock  profiles  traced  directly  from 
the  X-Y  recorder  plots.  The  maximum-slope  thickness  8  was  obtained  by  draw¬ 
ing,  by  eye,  a  straight  line  through  the  maximum  slope  point  of  the  profile  and 
measuring  the  distance  between  the  intersections  of  this  line  with  the  extrapo¬ 
lations  of  the  upstream  and  downstream  flow  regions  .  The  density  ratio  $2^1 
across  the  shock  was  determined  by  extrapolating  the  upstream  and  downstream 
densities  to  thp  midpoint  of  the  shock  It  is  believed  that  this  procedure 
will  account  properly  for  any  attenuation  of  the  beam  By  moving  the  shock- 
holder  choker  downstream  so  that  the  shock  was  swallowed,  measurement  of  the 
radiation  intensity  could  be  made  in  the  free  stream  and  compared  with  the 
known  density.  Some  attenuation,  of  the  order  of  15%  over  the  region  traversed., 
was  found.-  The  attenuation  depended  on  the  gas  density,  being  larger  at  higher 
gas  densities.  It  was  also  found  to  depend  on  the  alignment  of  the  spectrometer 
slit  image  with  the  beam,  and  this  would  vary  from  time  to  time-  After  some 
initial  attempts,  it  was  decided  not  to  correct  the  measurements  for  beam  atten¬ 
uation,  principally  because  such  corrections  should  not  make  any  measurable 
differences  in  the  maximum  slope  thicknesses,  density  ratios,  or  normalized 
shock  profiles. 

In  the  free  jet  flow,  there  is  a  small  negative  density  gradient  in 
the  free  stream  ahead  of  the  shock  waves,  although  it  is  hardly  discernable  in 
Figs,  3b  and  3d,  A  plane  shock  wave  in  a  divergent  stream  is  to  first  approx¬ 
imation  equivalent  to  a  curved  shock  wave  in  a  parallel  flow,  and  consequently 
in  the  free  jet  a  small  positive  density  gradient  is  produced  behind  a  shock 
wave-  The  effect  of  shock  curvature  itself  would  be  to  add  to  this  gradient, 
so  that  the  total  gradient  observed  is  due  to  a  combination  of  flow  divergence 
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and  shock  curvature,  plus  any  variation  in  radiation  intensity  associated  with 
attenuation  of  the  primary  beam  and  halo.  The  argon  trace  shown  in  Fig.  3c  is 
exhibited  as  an  extreme  case  of  downstream  density  gradient.  Most  argon  traces 
did  not  exhibit  this  much  downstream  gradient,  but  in  general  the  gradients 
behind  the  argon  shocks  were  greater  than  those  behind  the  nitrogen  shocks, 
possibly  due  to  the  greater  halo  intensity  surrounding  the  primary  beam  in 
argon.  Gradient  effects  were  much  less  behind  the  nozzle -produced  shock 
waves,  indicating  that  these  shock  waves  were  very  nearly  plane. 

A  very  curious  anomaly  was  observed  in  traces  obtained  for  the  free- 
jet  helium  shock  waves,  a  typical  one  of  which  is  shown  in  Fig.  3d.  In  all 
of  the  free- jet  measurements,  a  "dip"  in  radiation  intensity  was  found  in  the 
upstream  tail  of  the  shock  wave.  Whether  or  not  this  dip  corresponds  to  a  dip 
in  fluid  density  is  unknown.  It  was  absent  from  the  helium  free-jet  flow  when 
the  shock  was  swallowed  by  the  shock-holder,  and  it  was  never  observed  in 
argon  or  nitrogen.  It  was  likewise  not  observed  in  the  helium  shock  waves 
produced  in  nozzle  flows.  The  helium  free-jet  data  for  maximum-slope  thick¬ 
ness  are,  because  of  this  effect,  subject  to  greater  uncertainty  than  the 
argon  and  nitrogen  data,  and  the  situation  is  even  less  clear  with  regard  to 
experimental  density  ratios  for  the  helium  free-jet  data. 

Reciprocal  maximum  slope  thicknesses  for  argon  and  helium  are  plotted 


in  Fig.  4,  and  for  t>itrv/^cn  ir  F^-.  5.  Rather  than  use  the  upstream  mean  free 

path  we  have  chosen  to  normalize  5^  with  the  length  L  =  p  /pu  (thus 

*  •  *. 

forming  a  shock  thickness  Reynolds  number  Re^  =  pu  o  /|i  ).  There  are  sev- 

* 

eral  reasons  for  this  choice.  The  length  L  ,  being  based  on  the  viscosity 

* 

at  temperature  T  ,  is  much  more  representative  than  of  the  conditions 

in  the  middle  of  the  shock  wave  where  6  is  measured.  Another  advantage 

m 

which  accrues  from  the  use  of  L  is  that  theoretical  predictions  show  that 


is  quite  insensitive  to  the  viscosity- temperature  law.  and  also  that 
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L  /B  asymptotes  to  a  constant  value  as  M^  -»  «.  Finally,  and  most  important, 

were  we  to  use  as  a  reference  length,  we  would  be  required  to  evaluate 

the  viscosity  at  very  low  temperatures,  where  it  is  not  at  all  well  known. 

In  Fig.  4,  the  free  jet  data  for  all  stream  densities  at  a  given  Mach  number 

have  been  averaged  and  plotted  as  a  single  point.  Error  flags  of  ±7%  would 

encompass  all  of  the  data  at  each  Mach  number.  Individual  points  have  bean 

plotted  for  each  stream  density  for  the  nozzle  data,  because  the  Mach  number 

changed  with  stream  density  due  to  boundary  layer  effects  as  noted  previously. 

No  systematic  trend  of  with  was  noted  for  the  free  jet  experiments, 

indicating  that  in  the  range  of  densities  covered  by  these  experiments  the 

response  of  the  electron  beam  was  linear  with  density.  However,  there 

appeared  to  be  a  trend  of  increasing  6  with  p,  in  the  Mach  3  nozzle 

m  i 

data.  As  discussed  later  on,  this  may  be  associated  with  nonlinearity  of  the 
fluorescence-density  relationship. 

Figure  4  also  shows  a  comparison  between  our  data  and  the  predictions 

of  the  Navier-Stokes  and  Mott-Smith  theories.  The  Navier-Stokes  curve  was  taken 

20 

from  the  paper  of  Schwartz  and  Hornig,  while  the  Mott-Smith  curve  is  based  on 

21 

an  average  of  the  calculation  of  Muckenfuss  and  of  Schwartz  and  Hornig.  The 

22 

viscosities  for  the  experimental  points  were  calculated  from  the  Bromley-Wilke 
cables.  It  can  be  seen  that  our  data  agree  quite  well  with  the  Mott-Smith 
theory,  at  the  higher  Mach  numbers.  At  the  lowest  Mach  numbers,  our  data  fall 
between  the  Mott-Smith  and  Navier-Stokes  theories,  but  a  little  closer  to  the 
former.  Thus,  the  present  low  Mach  number  data  represent  shock  thicknesses 
a  little  greater  than  those  obtained  previously  by  Talbot  and  Sherman. 

For  purposes  of  comparison  we  have  in  Fig.  4  plotted  mean  curves 
representing  the  argon  data  obtained  by  Russell  and  by  Camac,  as  transposed 
to  L  /Bn  coordinates.  The  relationship  between  and  L  ,  for  Russell's 

data,  is 
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while  for  Camac's  data  it  is 


Both  formulas  give  essentially  the  same  result  for  7  =  5/ 3 o  Since  both  of 

these  investigators  worked  in  shock  tubes,  with  300°KS  the  corresponding 

values  of  T  at  large  Mach  numbers  go  as  high  as  5000°K,  where  viscosity  data 

are  subject  to  some  uncertainty .  For  the  high  tempera tures,  we  have  used  the 

23 

argon  viscosity  values  recommended  by  Amdur  and  Mason.  These  values  agree 
with  those  of  Bromley  and  Wilke  for  T  <  1200°K.>  but  above  this  temperature 
the  Amdur-Mason  viscosities  are  substantially  larger  than  the  Bromley-Wilke 
values. 

It  can  be  seen  that  Russell’s  data  and  our  own  are  in  reasonably 

good  agreements  while  Camac’s  experiments  yielded  slightly  thicker  shock  waves. 

Of  course,  a  substantially  different  viscosity-temperature  law  could  change 

the  relative  positions  of  the  data.  Although  we  have  not  shown  the  argon 

data  obtained  by  Linzer  and  Hornig  for  the  range  1.7  <  M^  <  5,  it  may  be 

noted  that  these  also  agree  quite  well  with  Mott-Smith  theory.  It  is  also 

of  interest  to  note  that  Russell's  data  were  in  quite  good  agreement  with 

24 

shock  wave  thicknesses  calculated  by  Chahine  and  Narasimha  and  Anderson 
25 

and  Macomber  on  the  basis  of  the  Bhatnagar-Gross-Krook  model  of  the 
Boltzmann  equation,  after  an  ad  hoc  correction  was  applied  to  change  from 
Prandtl  No.  =  1  to  Prandtl  No.  =  2/3.  Apparently,  if  this  Prandtl  Number 
correction  can  be  accepted,  there  is  little  difference  between  the  predic¬ 


tions  of  the  B-G-K  and  the  Mott-Smith  theories. 
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Reciprocal  shock  thicknesses  obtained  in  nitrogen  are  plotted  in  Fig, 

5,  following  the  same  scheme  as  that  used  for  the  monatomic  gas  data.  Also 
shown  for  comparison  are  the  data  of  1 Inzer  and  Hornig  and  Camac,  For  Camac's 
data,  T  extends  to  values  in  excess  of  400Q°R,  and  we  have  evaluated  ji(T  ) 
by  both  the  Sutherland  and  the  Bromley -Wilke  formulas.  The  latter  predicts  some¬ 
what  larger  viscosities  for  T  >  10003K,  It  can  be  seen  that  our  results  fall 
about  where  one  might  estimate  an  extrapolation  of  the  Linzer-Hornig  data  to  go. 
However,  our  measurements  yielded  shock  waves  considerably  thinner  than  those 
observed  by  Camac.  The  lack  of  agreement  between  any  of  the  data  and  Navier- 
Stokes  theory  with  bulk  viscosity  K  -  2^/3  is  evident.  The  Navier-Stokes 

curve  was  taken  from  Schwartz  and  Hornig;  K  ~  2ja/3  is  the  value  suggested  by 
2 

Sherman^s  low  Mach  number  measurements  in  air.  The  reason  for  the  rather 

large  discrepancy  between  Camac  s  data  and  our  own  is  not  clear.  It  would  be 

•x 

difficult  to  assign  the  ertire  discrepancy  to  the  uncertainty  in  |i(T  )  at 
high  temperature,  although  this  could  be  a  contributory  factor, 

3 . 3  Shock  Wave  Density  Ratios 

A  comparison  between  the  experimental  and  theoretical  density  ratios 

in  the  form  A  =  (p  /p  )  /(p^/p.)  provides  another  test  of  the  accuracy 

ti  x  sxp t  /  l  « neor 

of  the  electron-beam  fluorescence  technique.  In  the  Mach  1,5  nozzle,  both  the 
argon  and  nitrogen  data  corresponded  to  values  in  the  range  A  =  1  ±  0.02, 
whereas  the  rather  small  amount  of  helium  data  obtained  gave  A  =  0,94  ±  0.03, 

The  agreement  was  poorer  in  the  Mach  3  nozzle.  Here  the  nitrogen  data  gave  an 
average  A  =  0,95,  while  the  helium  data  gave  A  =  0,90,  The  argon  data,  on 
the  other  hand,  yielded  A  =  1,07,  We  believe  chis  poorer  agreement  in  the 
Mach  3  nozzle  is  due  principally  to  the  failure  of  the  linear  relationship 
between  beam  intensity  and  gas  density,  and  perhaps  complicated  by  halo  effects. 
The  absolute  free-strearo  gas  density  in  the  Mach  3  nozzle,  at  even  the  lowest 
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stagnation  pressures  in  both  argon  and  nitrogen,  barely  overlapped  the  highest 
free  stream  densities  in  the  Mach  1.5  nozzle  or  free  jet  tests.  If  this  ex¬ 
planation  is  correct,  we  can  infer  a  rough  upper  bound  estimate  for  the  maximum 
density  at  which  the  fluorescence  method  can  be  used  with  confidence,  at  least 
for  and  argon.  If  we  take  conditions  behind  the  shock  wave  as  representa¬ 
tive,  then  from  both  the  nitrogen  and  the  argon  data  this  upper  bound  on 

*7  3 

density  appears  to  be  in  the  vicinity  of  5  X  10  gm/cm  .  In  round  numbers, 
this  density  corresponds  to  pressures  of  about  0.25  Torr  in  argon,  and  0.35 
Torr  in  nitrogen  at  T  =  300°K.  Although  the  estimates  are  rough,  they  are 
not  in  conflict  with  Gadamer's  data. 

The  halo  radiation,  especially  in  argon,  may  also  contribute  to  the 
departures  from  A  =  1  in  the  Mach  3  nozzle.  This  effect  is  hard  to  estimate, 
because  we  do  not  have  any  reliable  measurements  on  the  intensity  and  extent 
of  the  halo  radiation  as  a  function  of  gas  density.  We  did,  however,  make 
several  jeans  of  shock  waves  with  the  spectograph  slit  imaged  so  as  to  observe 
only  the  halo  radiation  above  the  primary  beam.  These  shock  profiles  looked 
quite  similar  to  those  taken  along  the  primary  beam,  but  were  slightly  (about 
10%)  thicker.  Perhaps  this  accounts  for  the  fact  that  the  Mach  3  data  seem 
to  be  consistently  a  little  farther  below  the  Mott-Smith  curve  than  the  other 
data. 

The  experimental  density  ratios  for  the  shock  waves  produced  in  the 

free  jet  flows  can  also  be  compared  with  the  theoretical  values  of  the  modified 

Rankine-Hugoniot  conditions  for  low-Reynolds -number  spherical  shock  waves.  The 

low-Reynolds-number  correction  to  O^A^theor  due  to  shock  curvature  was 

26 

first  examined  by  Germain  and  Giuraud.  Numerical  results  for  a  spherical 

27 

shock  wave  have  been  calculated  by  Oberai.  Along  the  axis  of  symmetry  (see 
Fig.  6),  the  corrected  density  ratio  ^2^1^curv  :''s  accor<^nS  to  Oberai 

related  to  the  Rankine-Hugoniot  value  for  plane  shocks  (^^l^RH5  ^ 


*) 
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(p2/p!> 


=  (p 


curv- 


^RH  ( 


1  + 


P11  +  P22 


(3) 


The  factor  (p^,  +  p22^/p0  was  found  by  Oberai  to  have  the  values  -7.045  and 
-7.70  for  7  **  1.4  and  1.67,  respectively.  The  parameter  €  is  an  inverse 
Reynolds  number,  defined  by  6  =  p  /pu  R^,  where  R^  is  the  nose  radius  of 


the  body  producing  the  shock  wave,  but  we  have  ignored  the  small  difference 

*  / 

between  R^  and  the  shock  wave  radius  Rg  and  have  used  instead  €  =  p  /pu  Rg 
The  theory  makes  use  of  hypersonic  approximations,  and  is  therefore  only 


applicable  for,  say,  M  ^  3. 

For  a  curves  shock  wave  in  a  slightly  divergent  source-flow,  the 
effective  shock  radius  one  must  use  in  evaluating  e  is  to  first  approxi¬ 
mation 


R 
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where  R  is  the  radial  distance  from  the  shock  to  the  source,  as  shown  in 
o 

Fig.  6.  We  assume,  however,  that  our  shock  waves  are  plane,  so  we  have  taken 

R  =  R  .  In  Fig.  7  we  have  plotted  values  of  A  =  (p9/p..)  /( p9/p.)  and 

s  °  z  1  exp  1  1  RH 

)  for  our  free-jet  data.  One  observes  that 

curv. 

when  the  curvature  correction  is  included,  all  the  data  can  be  represented 
hy  ^corr  =  0*99  ±0.02,  whereas  without  the  correction  A  =  0.97  ±0.02. 

(We  have  corrected  the  low  Mach  number  data  as  well,  even  though  the  theory 
is  a  hypersonic  one.)  Evidently,  the  curvature  correction  improves  the  agree¬ 
ment  with  theory  to  a  point  where  the  values  of  A  obtained  in  the  free-jet 
experiments  agree  very  well  with  those  obtained  in  the  Mach  1.5  nozzle,  and 
one  may  perhaps  regard  this  as  a  tentative  verification  of  the  theory. 

Clearly,  a  more  detailed  investigation  of  this  problem  would  be  desirable 
and  is  within  the  capability  of  the  electron-beam  fluorescence  method. 


corr 
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The  density  gradient  behind  the  shock  wave,  due  to  flow  divergence  and 
shock  curvature,  can  easily  be  calculated  from  elementary  inviscid  flow  theory, 
if  the  shock  wave  is  assumed  to  be  vanishingly  thin.  One  finds 


^2  2M22p2  r  IQi*  •  1)  1  (l_  v 

“  1  -  M22  L  2  *  (r-DMj2  j  '  v  \  \) 


(5) 


We  compared  this  prediction  with  the  experimental  density  gradients,  assuming 
Rg  =  0® «  but  the  result  was  inconclusive.  The  experimental  values  were  of  the 
same  order  as  the  theoretical  ones,  but  the  scatter  was  as  much  as  100%.  The 
scatter  could  be  due  to  a  variety  of  effects,  such  as  noise  in  the  traces, 
small  deviations  from  planarity  of  some  of  the  shock  waves,  and  halo  and 
scattering  effects. 


4.0  CONCLUSIONS 

Shock  wave  thicknesses  and  density  ratios  were  measured  in  helium,  argon 
and  nitrogen  by  means  of  the  electron  beam  fluorescence  method,  over  the  Mach 
number  range  1.5  <M<-  17.4,  in  a  low  density  wind  tunnel.  Both  contoured-nozzle 
and  free-jet  flows  were  used.  The  shock  thicknesses  in  argon  and  helium  agreed 
quite  well  with  Moct-Smith  theory  at  the  higher  Mach  numbers,  and  were  between 
Navier-Stokes  and  Mott-Smith  theory  at  the  lowest  Mach  number.  For  the  higher 

Mach  numbers  good  agreement  was  obtained  with  other  measurements  made  by  the 

7  9  1  20 
electron  beam  attenuation,  electron  beam  scattering,  and  optical  reflectivity  8 

methods.  It  is  concluded  that  the  electron  beam  fluoresence  technique  is  at 

least  as  accurate  as  any  of  these  other  methods.  In  nitrogen,  the  measured 

shock  thicknesses  were  considerably  greater  than  the  predictions  of  the  Navier- 

Stokes  theory.  They  agreed  well  with  previous  data  obtained  by  the  optical 

reflectivity  method,  but  were  considerably  thinner  than  these  measured  by  the 

electron  beam  scattering  technique. 
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Measured  density  ratios  across  the  shock  waves  in  nitrogen  and  argon  were 
in  very  good  agreement  with  theory  in  a  nominal  Mach  1.5  nozzle.  In  a  nominal 
Mach  3  nozzle,  agreement  was  poorer,  the  reason  apparently  being  a  failure  of 
the  linear  relationship  between  radiation  intensity  and  gas  density.  This 
interpretation  leads  to  rough  upper  bounds  of  p  =  0.25  Torr  and  p  =  0.35  Torr 
at  room  temperature  for  the  pressures  in  argon  and  nitrogen,  respectively,  at 
which  the  fluorescence  method  can  be  used  with  confidence. 

The  density  ratios  across  the  shock  waves  in  the  divergent  free-jet  flows 
were  subject  to  a  viscous  correction  due  to  shock  wave  curvature.  When  a 
theoretical  correction  was  applied,  agreement  between  experiment  and  theory 
was  very  good,  which  can  be  interpreted  as  a  tentative  verification  of  the 
theory . 

An  anomalous  dip  was  observed  in  radiation  ahead  of  the  helium  shock 
waves  produced  in  the  free-jet  flows.  This  dip  was  not  observed  in  the 
helium  shock,  waves  produced  in  the  nozzles,  nor  in  any  of  the  argon  or 
nitrogen  shock  waves.  The  reason  for  the  dip  is  not  understood. 
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FIG.  1.  SCHEMATIC  OF  EQUIPMENT  FOR  SHOCK  WAVE  MEASUREMENTS  IN  FREE  JET 
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FIG.  2b.  BEAM  RADIATION  PROFILES  IN  HELIUM 


Mach  4.17,  30  kV  beam  at  2.0  mA 

Static  density  equivalent  0.35  Torr  at  297'K 

(Note  scale  difference  for  the  3*P  and  4^)  lines) 
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FIG.  3b.  SHOCK  PROFILE  IN  NITROGEN.  FREE  JET,  Mj  -  9.45,  Pn  “  42 
3  X  10'8  A  full  scale.  Beam  30  kV  at  1.5  tnA. 


Torr, 
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FIG.  3c.  SHOCK  PROFILE  IN  ARGON.  FREE  JET,  Mj  =  17.4,  pQ  =  61.2  Torr, 


3  X  10  A  full  scale.  Beam  30  kV  at  2.5  raA. 
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FIG.  3d.  SHOCK  PROFILE  IN  HELIUM.  FREE  JET,  Mj  -  17.4,  pQ  =  51  Torr, 
3  X  10‘9  A  full  scale.  Beam  30  kV  at  4.0  mA. 
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FIG.  4.  MAXIMUM-SLOPE  THICKNESSES  FOR  ARGON  AND  HELIUM 
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FIG.  5.  MAXIMUM  SLOPE  THICKNESSES  FOR  NITROGEN 
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FIG.  6.  NOMENCLATURE  FOR  SHOCK  WAVE  CURVATURE 
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FIG.  7.  FREE  JET  SHOCK  HAVE  DENSITY  RATIOS  WITH  AND  WITHOUT  CURVATURE  CORRECTIONS 
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